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The description of realistic strongly correlated systems has recently advanced through the combi- 
nation of density functional theory in the local density approximation (LDA) and dynamical mean 
field theory (DMFT). This LDA+DMFT method is able to treat both strongly correlated insulators 
and metals. Several interfaces between LDA and DMFT have been used, such as (N-th order) Linear 
MufSn Tin Orbitals or Maximally localized Wannier Functions. Such schemes are however either 
complex in use or additional simplifications are often performed (i.e., the atomic sphere approxi- 
mation). We present an alternative implementation of LDA-I-DMFT, which keeps the precision of 
the Wannier implementation, but which is lighter. It relies on the projection of localized orbitals 
onto a restricted set of Kohn-Sham states to define the correlated subspace. The method is imple- 
mented within the Projector Augmented Wave (PAW) and within the Mixed Basis Pseudopotential 
(MBPP) frameworks. This opens the way to electronic structure calculations within LDA-I-DMFT 
for more complex structures with the precision of an all-electron method. We present an application 
to two correlated systems, namely SrVOa and /3-NiS (a charge-transfer material), including ligand 
states in the basis-set. The results are compared to calculations done with Maximally Localized 
Wannier functions, and the physical features appearing in the orbitally resolved spectral functions 
are discussed. 
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I. INTRODUCTION 



The description of strong electronic correlations in a 
realistic framework has become an issue of major im- 
portance in current condensed matter research. Due to 
the fast progress in the preparation of novel materials, 
especially in effectively reduced dimensions, and the ad- 
vances in experimental techniques in order to probe such 
systems, providing an adequate theoretical formalism 
that can handle explicit many-body effects in a material- 
specific context is crucial. It is evident that standard 
density functional theory (DFT) in the local density ap- 
proximation (LDA) cannot meet those demands and, at 
least, has to ally with manifest many-body techniques. 
In this respect, the combination of DFT-LDA with the 
Dynamical Mean-Field Theory (DMFT), the so called 
LDA-I-DMFT approach, has proven to be an invaluable 
method to face the challenge. Already numerous stud- 
ies within the LDA-I-DMFT framework have shown that 
this theory is capable of describing the effects of strong 
correlations in a realistic context, such as Mott transi- 
tions and volume-collapse transitions in d and /-electron 
systems, effective mass enhancement, local moment for- 
mation and magnetism, spectral-weight transfers, orbital 
physics, etc... 

In view of these successes, it is crucial to push further 
the range of applicability of LDA-I-DMFT methods. To 
this aim, implementing LDA-I-DMFT within those elec- 
tronic structure methods which are highly accurate and 
allow for the possible treatment of larger systems, is cer- 
tainly an important endeavor. One of the main aims 



of this paper is to report on such new implementations, 
within a projector augmented-wave (PAW) and a mixed- 
basis pseudopotential (MBPP) method. The formalism 
used in this paper is actually quite general and allows for 
an implementation of LDA-f DMFT in a very large class 
of electronic structure methods. We present this formal- 
ism in an arbitrary basis set, which should make this task 
easier for other implementations. However, in our specific 
implementation, we make an extensive use of the Bloch 
basis set of Kohn-Sham orbitals. In order to implement 
LDA-I-DMFT, we use local orbitals constructed by pro- 
jecting atomic-like orbitals onto a restricted set of Bloch 
states, a strategy similar to that introduced by Anisimov 
et alX. 

Another important route is to investigate how to op- 
timize the application of LDA+DMFT for a specific ma- 
terial, or class of materials of interest. Indeed, it is im- 
portant to realize that, while the results of course do 
not depend on the specific basis set used in the calcu- 
lation, they will actually depend on the specific set of 
local orbitals for which many-body effects will be taken 
into account, or more precisely (following the terminol- 
ogy introduced in 0]) on the "correlated" subspace C of 
the full Hilbert space spanned by these orbitals. Indeed, 
the DMFT treatment applies local interaction terms to 
those orbitals only, and furthermore neglects all non-local 
components of the self-energy. This notion of locality is 
defined with respect to the specific choice of the local or- 
bitals defining C, and the accuracy of the DMFT approx- 
imation cannot be expected to be identical for different 
choices. 
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Early implementationgMi^ of LDA+DMFT used the 
linear muffin-tin orbital^ (LMTO) framework, and the 
correlated orbitals were frequently identified with a spe- 
cific subset of the LMTO basis functions, having d- or 
/-character. There is no specific reason to pick the local 
orbitals as members of the basis set however, and a set 
of atomic-like functions may prove to be better from the 
physical point of view. Recently, several works have used 
different kinds of Wannier functions (WFs) as correlated 
orbitals, starting with the work of Pavarini et al^ us- 
ing Nth-order muffin-tin orbitals- (NMTO). Anisimov et 
al} used WFs constructed from a projection on a subset 
of the Bloch functions^ii^, and Lechermann et al. used 
maximally-localizedii WFs and compared the results to 
those using NMTOs. 

The energy window spanned by the basis functions 
which are retained in the implementation, and the spa- 
tial extension of the local orbitals defining C are impor- 
tant physical issues for the description of a given ma- 
terials. Those issues become particularly important for 
charge-transfer materials (e.g. late transition-metal ox- 
ides, sulfides or selenides) in which the ligand states must 
be kept in the basis set in order to reach a satisfactory 
physical picture. In the present article, we address these 
issues and provide explicit comparisons between calcula- 
tions performed with more spatially extended local or- 
bitals (hence spanning a smaller energy window when 
projected onto Bloch functions) and more localized local 
orbitals (hence spanning a wider energy window) . 

This article is organized as follows. Sec. |ll]is devoted 
to a presentation of the general theoretical framework. 
Sees, mil and IIVI are devoted to applications to two com- 
pounds, which are considered as tests for the method: the 
transition- metal oxide SrVOa in Sec. IIIII and the charge- 
transfer sulfide NiS in Sec. IIVI 



II. THEORETICAL FRAMEWORK 

A. LDA+DMFT formalism in an arbitrary basis 

set 

In order to implement DMFT within realistic elec- 
tronic structure calculations of correlated-electron ma- 
terials, one has to set up a formalism which keeps 
track of the real-space (i.e., quantum chemical) and the 
reciprocal-space (i.e., solid-state) aspects on an equal 
footing, while being computationally efficient. Follow- 
ing the work of Lechermann et al^, we therefore distin- 
guish between the complete basis set {liJka}} in which 
the full electronic structure problem on a lattice is for- 
mulated (and accordingly the lattice Green's function is 
represented), and local orbitals \Xm) which span a "cor- 
related" subspace C of the total Hilbert space. Many- 
body corrections beyond LDA will be considered inside 
this subspace. The index a labels the basis functions for 
each wave vector k in the Brillouin zone (BZ). The in- 
dex R denotes the correlated atom within the primitive 



unit cell, around which the local orbital is centered, 
and m=l, . . . , M is an orbital index within the correlated 
subset. Projection on that subset, for atom type R, is 
done with the projection operator: 
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In Ref. it was shown that the DMFT self-consistency 
condition, which relates the impurity Green's function 
G-Sip to the Green's function of the solid computed lo- 
cally on atom R, reads: 
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denotes the Bloch transform of the local orbitals whereby 
T denotes the Bravais lattice translation vectors. Note 
that the object in the second line of Eq. ^ is, of 
course, nothing else than the full lattice Green's function 
Gqcj' (k, icjn) in the chosen {|-Bka)} basis. The Kohn- 
Sham (KS) hamiltonian HKs(k) is obtained by solving 
the Kohn-Sham equations, which yields eigenvalues eki/ 
and Bloch wave- functions I'I'kiy) {1^ is the band index). It 
can be expressed in the {l-Bka)} basis set as: 
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In order to obtain the self-energy for the full solid, one has 
to promote ("upfold") the DMFT impurity self-energy 
to the lattice via^: 
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whereby a double-counting correction S^„/ , taking care 
of correlation effects already accounted for in the LDA 
hamiltonian (see appendix [B|) has to be included in the 
general case. 

The above equations form the general LDA+DMFT 
framework, in a general arbitrary basis set. Any specific 
implementation must then make a choice for: 

i) The set of local orbitals {|Xm)} spanning the cor- 
related subspace 

ii) The specific basis set {|i?ka)} in which these equa- 
tions will be implemented 

It is important to realize that the results will certainly 
depend on the specific choice of {|Xm)}- the quality of 
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the DMFT approximation will indeed depend on how 
the local orbitals are picked such as to minimize non- 
local contributions. In contrast, for a given choice of 
{|Xm)}'s, the results should be in principle independent 
of the basis set {|i?ka)} which is chosen for the imple- 
mentation. However, in practice, considerations of nu- 
merical efficiency do limit the size of the basis which 
can be handled. Indeed, Eq. ^ involves inverting a 
matrix (at each k-point and for each frequency) of size 
Nh X Nh, in which TVf, is the number of basis functions 
which are eventually retained. Hence, in practice, one 
will restrict the basis set to a certain set W of bands, 
as will be described in more details below. Regarding 
the choice of local orbitals, we recall that, in Ref. 0, 
different kinds of Wannier functions were used and com- 
pared to one another, namely maximally-localized^^ and 
A^th-order linear- mufhn-tin orbitals (NMTO)^. The con- 
struction of such functions requires rather sophisticated 
procedures. In the present work we use a somewhat 
lighter implementation, with the same demands on accu- 
racy, by constructing the {|Xm)}'s out of entities that are 
already existing in most of the common band-structure 
codes, namely using the decomposition of local atomic- 
like orbitals onto the basis function retained in the set 
W. This is very similar in spirit to the construction pro- 
posed by Anisimov et al} in the LMTO framework. The 
present construction, and the basic equations developed 
in this section, are fairly general however and can be im- 
plemented in an arbitrary electronic structure code. 



B. LDA+DMFT formalism in the Bloch basis set 

In this section, we focus on a very simple choice for the 
basis set, namely the Bloch basis itself {|i?kQ)}={|^'ki/)}- 
This basis is most conveniently used, since it is a direct 
output of any DFT-LDA calculation and furthermore di- 
agonalizes the KS hamiltonian Hj^3(k) = S,^,^'e\^i,. 

The basic LDA+DMFT equations in the Bloch basis 
set are easily written up, using the projection matrix 
elements of the local orbitals onto the Bloch functions, 
defined as: 

n^.(k) ^ (xLl^k.) , P^;A^) ^ (*k.|xD • (6) 
Equations ©J® then read: 

Glt'^M = ^P„^(k)G^|„(k,z^„)F,^;,,(k) (7) 
AE^l„(k,zc.„) = ^^P,^^*(k)AE3,(z^„)^™v'P)) 

R, mm' 

where 

G",(iw„) = |[(iw„ + fJ-~ ekiy)^i/i/' - ASbKk, iw„)]"^ 
^^3' i^^n) = E™P, (iw„) - . 



C. Truncating the Bloch basis set, and choice of 
local orbitals 

As pointed out above, it is computationally impossi- 
ble to implement these LDA-I-DMFT equations without 
restricting oneself to a finite subspace of Bloch func- 
tions. Those Bloch functions span a certain energy win- 
dow, corresponding to a subspace W of the total Hilbert 
space. Naturally, the local atomic-like orbitals |Xkm) 
will in general have a decomposition involving all Bloch 
bands. The Bloch-transform of these local orbitals read: 

IXkJ-E(*k.|xL>l*k.) . (11) 

Note that, starting from an orthonormalized set of local 
orbitals (XmklXm'k') f^mm'feR'Jkk', it is easily checked 
that the matrix (^kiy|x^m) is unitary (from the complete- 
ness of the Bloch basis). Hence the Xm's can formally be 
viewed as Wannier functions associated with the com- 
plete basis set of all Bloch states. This property no longer 
holds, however, when the sum in pT|) is restricted to the 
subset W of Bloch band. Defining: 

IXkJ^ E<*k.|XkJI*k.) (12) 

it is seen that the functions |x^m) ^'^'^ not true Wannier 
functions associated with the subspace W since the trun- 
cated projection matrix is no longer unitary. However, 
these functions can be promoted to true Wannier func- 
tions Iw^) by orthonormalizing this set according to: 

K™)= E ^™n'(k)IXk™'> , (13) 

R'rn' 

where (k) is given by the inverse square root of the 
overlap matrix between the Wannier-like orbitals, i.e., 

O^^'(k) ^ (xLlx^™'> = E ^™^(k)Pj^™'(k)(14) 

5^^,(k) ^ |(o(k))-i/n (15) 

^ mm' 

Naturally, the functions are more extended in space 
than the original atomic-like functions Xm since they can 
be decomposed on a smaller number of Bloch functions, 
spanning a restricted energy range. 

In the end, LDA-f DMFT is implemented by taking for 
C the correlated subset generated by the set of functions 
jw^j). Since those functions have a vanishing overlap 
with all Bloch functions which do not belong to the set 
W, the LDA+DMFT equations can now be put in a com- 
putationally tractable form, involving only a NbX Nb ma- 
trix inversion within the selected space W. Hence, the 
equations which are finally implemented read: 

Gm^^i^^n) = E PmA^)Gllik,^u;^)PX'mQ) 
^p^l„(k,zc.„) = EE^-"(k)^S3(zc.„)P;^,,,(kT) 

R. mm' 
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where 

P^,(k) ^ ^w(k)P5.(k) , (18) 
F,^:(k) EE ^ 5;^^*(k)P5;,(k) . (19) 

R,'m' 

It is important to realize that the truncation to a Hmited 
set of Bloch functions was not reached by simply neglect- 
ing matrix elements between the local orbitals and Bloch 
functions outside this set, but rather by constructing a 
new set of (more extended) local orbitals such that the 
desired matrix elements automatically vanish, hence re- 
defining C accordingly. In this view, the choice of C and of 
W, although independent in principle, become actually 
inter-related. 

We also note that it is not compulsory to insist on form- 
ing true Wannier functions out of the (non-orthogonal) 
set Ix^m)- perfectly legitimate formally to choose 

the correlated subspace C as generated by orbitals hav- 
ing a decomposition in W, but not necessarily unitarily 
related to Bloch functions spanning W. Although or- 
thogonality of the Xm's is also not compulsory, several 
(but not all) impurity solvers used within DMFT do re- 
quire however that the x'm^ be orthogonal on a given 
atomic site. One possibility, for example, is to orthonor- 
malize this set on identical unit-cells only, i.e requiring 
that l^xm) ^^'^ l^T'm') space are orthogonal for 

T = T', but not in neighboring cells T ^ T'. This 
amounts to orthonormalize the set with respect to 

the k-summed overlap matrix, instead of the one com- 
puted at each k-point. 

In our actual implementations, the wave- functions 
spanning the correlated subspace C are obtained by fol- 
lowing the above orthonormalization procedure, starting 
from atomic-like orbitals Xm centered on the atomic site 
R in the primitive unit-cell. These local orbitals are ei- 
ther all-electron atomic partial waves in the PAW frame- 
work, or pseudo-atomic wave functions when using the 
MBPP code. Since, in the present work, we are not 
dealing with full charge self-consistency including self- 
energy effects, the matrix elements P^^Ck) and the wave- 
functions can be computed once and for all at the 
beginning of the DMFT cycle^. Details on the specific 
construction of the local orbitals used in this article, and 
the corresponding calculation of ([6|) are summarized in 
appendix |^ 

D. Physical considerations on the choice of the 
correlated subspace C and of the Wannier/Bloch 
space W 

Let us now discuss some physical considerations re- 
garding the choice of the truncated Bloch space W when 
using the LDA-I-DMFT framework to describe a given 
material. Operationally, this means that a certain num- 
ber of Bloch bands Nb (spanning a certain energy win- 
dow) will be retained when solving Eqs. (|16ll7p . As dis- 



cussed above, the choice of W also influences the actual 
definition of the correlated subspace C, since we require 
that the orbitals generating C can be expanded upon the 
basis functions generating W. 

Let us consider, to be specific, the case of a transition- 
metal oxide, such as for example SrVOa. This material, 
which is described in more details in the following section, 
has a set of three t2g bands, well separated from both 
the 02p and Cg bands, and containing nominally one d- 
electron. We can make two rather extreme choices when 
describing this material with LDA-I-DMFT: 

i) Focus only on a very limited set of low-energy Bloch 
bands, such as the three t2g bands, and generate W just 
from the three corresponding Bloch functions. In this 
case, we shall have also C — W, and the |wx„) will be 
Wannier functions unitarily related to these three Bloch 
bands. Since these bands span a narrow energy win- 
dow, this also means that these Wannier functions will 
be rather extended spatially: although centered on Vana- 
dium atoms, they will have a sizeable contribution on 
neighboring oxygen atoms as well. This kind of approach 
has been emphasized and studied in details in Ref. 0. 
Of course, it is then out of the scope of such approaches 
to investigate the indirect effects of correlations on bands 
other than the t2g ones. 

ii) Alternatively, one may choose a large energy win- 
dow to define W, including in particular all Bloch bands 
corresponding to 02p, V-t2g and Y-Cg. Then, the or- 
bitals Iwxm) defining the correlated subspace C may be 
chosen as having a component on Bloch states spanning 
a much larger energy range. As such, they will be more 
localized spatially, i.e closer to (vanadium) atomic-like 
orbitals. When working with such an enlarged space W, 
the physics of 02p and Cg states can also be addressed. 

One of the goals of the present paper is to present and 
compare calculations done with such different choices of 
W and C. Of course, in a fully first-principle approach, 
the screened interaction matrix elements should also be 
calculated (e.g in a GW framework) in a manner which 
is consistent with these choices of Hilbert spaces. This is 
left for future investigations however, and in the present 
work, these matrix elements will be taken as parameters. 

Let us note that the local orbitals are constructed in 
the present paper from an atomistic point of view. Hence 
there is no "entangling problem" as the one encountered 
when constructing maximally localized Wannier func- 
tions for strongly hybridized band complexes. Of course, 
in the present formalism, the Wannier functions obtained 
by projection of atomic orbitals onto Bloch states be- 
longing to a narrow energy window are not maximally- 
localized in the sense of Ref. [ll|, but this feature does not 
bring essential differences in the results, as clear from the 
results reported below. 

Finally, in order to relate the LDA-I-DMFT results to 
experiments performed using e.g, photoemission spec- 
troscopy, the real-frequency spectral functions must be 
obtained. This can be done using e.g. a maximum en- 
tropy treatment of Monte Carlo data, but is important 
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to understand how to connect the calculated quantities 
to physically observable spectra. The most direct out- 
put of the LDA+DMFT calculation is the local impurity 
spectral function, obtained (for orbital m) as: 

^ --ImG™P('^ + iO+) (20) 

TT 

This also corresponds to matrix elements of the full 
Green's function of the solid within atomic-like orbitals 
Xm's- This however, is not a quantity which can be eas- 
ily related to photoemission experiments, since the Xm's, 
when very localized spatially, extend over a large energy 
range. When considering a certain energy, contributions 
of other electrons with a different orbital character than 
the Xm's will contribute significantly to the photoemis- 
sion signal. Instead, if one is interested in the measured 
spectral function in a given energy window, one must 
consider the matrix elements of the full Green's function 
within Bloch (or Wannier) functions spanning that en- 
ergy range, namely: 

A.ico) = --ImV G^Uk,c^ + iO+) (21) 

k 

This quantity can be obtained either by first reconstruct- 
ing the local self-energy on the real-frequency axis by 
analytical continuation of the impurity Green's func- 
tion, or by direct analytic continuation of the imaginary- 
frequency Bloch Green's function G'''(k, 

III. APPLICATION: SRVO3 

SrVOs is a tlgCg metal. It is a good test case 
for LDA-fDMFT calculations because it is cubic and 
non magnetic and also the t2g bands are isolated from 
both Cg and oxygen p bands - in the LDA bandstruc- 
ture. Numerous calculations (including LDA+DMFT) 
and experiments have been done on this compound 
22i3^^a6i7a8i9^^^^^^^ Frojn tl^eggg gt^dies, it 

appears that there is a need to include correlations to 
describe correctly this compound. This is thus an ideal 
system to benchmark this new implementation. 

A. LDA 

For the PAW calculations, semicore states of V and 
Sr are treated in the valence. Valence states for Sr, V 
and O thus include respectively 4s4p5s, SsSpAsSd and 
2s2p states. PAW matching radius are 1.52 a.u., 1.92 
a.u. and 2.35 a.u. respectively. The experimental cubic 
crystal structure is used (space groupe Pm3m with lat- 
tice constant of 7.2605 a.u.). Atomic data are generated 
using ATOMPAW^i^. Calculations are done with the 
PAW code ABINIT^i29. The density of states and the 
LDA bandstructure are shown on Fig. [T] and [21 The pro- 
jection of the density of states on 0-p, \-t2g and V-e^ 
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FIG. 1: LDA band structure for SrVOa. 
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FIG. 2: (Color online) LDA total and projected density of 
states for SrVOa. 

and the character of the bands (see Fig. [3]) show that 
bands with Oxygen p and with Vanadium t2g characters 
are indeed isolated from the others. The hybridization 
between Oxygen and Vanadium orbitals is nevertheless 
clearly seen. 



B. LDA+DMFT 

In this section, we will present the results of our new 
scheme based on the projection of Bloch states upon Lo- 
cal Orbitals (Projection on Local Orbitals (PLO)). As 
previously emphasized, the number of KS bands used for 
the projection have to be chosen in a given range of en- 
ergy. The extension of the Wannicr-like renormalized or- 
bitals (Eq. [T5)l will depend on this choice. We will use the 
fact that the band structure of SrVOa is made of isolated 
blocs. The choices of W are summarized in table IIIIBI 
We will compare the results with LDA+DMFT calcula- 
tions done:^ with Maximally Localized Wannier Functions 
(MLWF) (see appendix [C]) . 
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FIG. 3: (Color online) LDA band structure for SrVOs com- 
puted in PAW, with "Fatbands" to show the amplitude of the 
projection of each band on a given atomic orbital (0-p, V-t2g 
and V-Bg). 



A'";, Nimp Energy range (eV) 
3 3 -1.5^1.8 
12 3 -8.0^1.8 
21 5 -8.0^13.0 



i;: 




FIG. 4: (Color online) Local impurity spectral function of 
SrVOs for U=4 eV within LDA+DMFT using the MLWF 
basis and the PLO scheme {Ni, — 3). Inset: Green's function 
in imaginary time. Note that when using such a small {t2g) 
energy window, the impurity spectral function and Bloch- 
resolved spectral function coincide. 

the same physical content and that PLOs constructed 
from a small energy window give very similar results to 
maximally localized Wannier functions^. 



TABLE I: (Color online) Different choices for the number of 
KS states Nt used within the PLO scheme of LDA+DMFT. 
The energy range spanned by these bands is also indicated. 
Nb — 3 corresponds to only the three bands of mainly 129 
character. Ni, = 12 corresponds to the 9 Oxygen-p bands 
and the 3 Vanadium-t29 bands (see Fig. |3}. Nt = 21 corre- 
sponds to a large number of bands, including eg bands and 7 
bands above. Aimp is the number of orbitals in the impurity 
model. When used, MLWF will be extracted with the same 
parameters. 



Calculations have been carried out at T=0.1 cV. A 
density-density interaction vertex is used'^°, with J=0.65 
eV, similarly to another study^. The Hubbard parame- 
ters U=4 eV and U=6 eV have been used in the calcu- 
lations. The impurity problem of DMFT is solved with 
Hirsch-Fye Quantum Monte Carlo with 128 time slices. 
The Around Mean Field formulation of the double count- 
ing is used in these calculations (see appendix |B| . 



1. t2g bands in the basis for W (Nb ~ 3). 

Spectral functions obtained from the PLO and MLWF 
schemes are reproduced on figure [H 

The impurity Green's function obtained within the two 
schemes are in very good agreement: they both show a 
lower Hubbard band at -1.8 eV and a upper Hubbard 
band at 2.5 eV. This shows that the two schemes contain 



2. Oxygen-p and Vanadium-t2g bands in the basis for W 
(Nb = 12;. 

Oxygen p states are now included in the calculation. 
MLWFs of t2g symmetry will thus be more localized, and 
will have less 0-p character, because the MLWFs of Oxy- 
gen p symmetry will have mainly the weight on Oxygens 
atoms. In the PLO scheme also, the Wannier orbitals will 
include more Bloch states, and thus will be closer to to 
localized local orbitals Xm- In other words, because the 
Wannier function of the local impurity problem is now 
constructed from a wider energy windows which includes 
oxygen p states, the corresponding spectral function will 
have a non negligible weight in the energy area of the 
oxygen-p states. In the LDA case without Hubbard cor- 
rection, this spectral function would in fact be identical 
to the projected t2g density of states plotted on Fig[T] 

In order first to benchmark our calculation using PLO 
with respect to calculations using MLWF, we plot on fig- 
ure [5] the comparison of the impurity Green's function 
in the two cases. The local impurity spectral function 
shows a large hybridization band in the area of the oxy- 
gen p states. We emphasize that this band is mainly not 
a many-body feature: it is a manifestation of hybridiza- 
tion with oxygen states, and is readily visible at the pure 
LDA level. The lower Hubbard band is in fact partly con- 
tained in this spectral function as a shoulder in the hy- 
bridization band around -1.5 eV. However, disentangling 
the Hubbard band from the hybridization contribution 
is somewhat difficult to achieve with the Maximum en- 
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U=4 eV, J=0.65 eV 



FIG. 5: (Color online) Local impurity spectral function of 
SrVOa for U=6 eV within LDA+DMFT using the PLO 
scheme (A^;, = 12). Inset: Green's function in imaginary time 
for both the PLO and the MLWF schemes. 



tropy method. The upper Hubbard band is visible in the 
impurity spectral function, as a hump around 3.0 eV. 

These many-body features (lower and upper Hubbard 
band) are more clearly revealed, however, by looking at 
the (k-averaged) spectral function of t2g Bloch states 
PT|) . It is also useful to look at the spectral functions 
of the other Bloch states in the basis-set. The summa- 
tion of these spectral functions over bands belonging to 
the same group (e.g Oxygen-p or Vanadium-t2g states, 
ie bands with mostly Oxygen-p or Vanadium-t2£, char- 
acter ) enables us to have a clear view on the impact 
of correlation on LDA bands. The spectral function of 
Oxygen-p states, and Vanadium-t2g states are plotted 
on figure [51 We emphasize that these spectral functions 
do not correspond to atomic-like orbitals with Oxygen-p 
character and Vanadium-t2g character (the latter being 
the impurity spectral function plotted on Figure O. 

The Hubbard band appears as a hump in the k- 
averaged Bloch spectral function corresponding to the 
t2g states for U=4 eV. This hump is located between the 
0-p states and the t2g-states quasiparticle peak. The 
Hubbard band is more clearly resolved for U=6 eV. In 
this case, however, it is hidden inside the Oxygen p band. 
We shall come back to this point at the end of this sec- 
tion. 

The fact that a higher value of U is necessary in this 
case (with respect to iVf, = 3.) to recover the lower Hub- 
bard band is consistent with the fact that Wannier func- 
tions are more localized for Nt — 12. 



3. Large number of bands in the basis for W (Nt — 21) 



In this case, the Cg states are included in the calcu- 
lation. The impurity model is now solved with all five 
d— orbitals. The agreement between impurity Green's 
functions computed in the MLWF and PLO schemes is 
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FIG. 6: (Color online) Total Bloch spectral function for 
SrVOa, and spectral functions of Oxygen-p, and Vanadium- 
t2g states for \]—4 eV (a) and U=6 eV (b) within 
LDA-I-DMFT using the PLO scheme (iVb = 12). These spec- 
tral functions are not the local orbitals projected spectral 
functions (see text): the t2g local impurity spectral functions 
are plotted on figure (5] 



shown on figure [T] Note again, that resolving the Hub- 
bard band from the impurity Green's function is diffi- 
cult because this quantity is dominated by hybridization 
effects with oxygen states. Again, we have to turn to 
Bloch-resolved spectral functions, plotted in Fig. [5] for 
Oxygen-p, Vanadium-t2g and Vanadium-Cg Bloch states. 
The results are quite similar to the previous ones with 
Nf, = 12, with a lower Hubbard band clearly visible for 
the t2g states at about -2.0 eV in the spectral function. 
Note that here we use U — 6 eV. It shows that the basis 
of Kohn-Sham bands is adapted to the calculation: the 
convergence of physical properties as a function of the 
number of bands is rather fast. 
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FIG. 7: (Color online) Local impurity spectral function of 
SrVOa for U=6 eV within LDA+DMFT using the PLO 
scheme {Nf, — 21). Inset: Green's function in imaginary time 
for both the PLO and the MLWF schemes. 
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FIG. 8: (Color online) Total Bloch spectral function for 
SrVOa, and spectral functions of Oxygen-p, Vanadium-t2g 
and Vanadium-eg states for U=6 eV within LDA+DMFT us- 
ing the PLO scheme (A^i, = 21). These spectral functions are 
not the local orbitals projected spectral functions (see text): 
the t2g and eg local impurity spectral functions are plotted 
on figure [T] 



4- Discussion 

The comparisons made above should not be simply 
thought of as a convergence study as a function of the 
size of the Bloch basis W. Indeed, as Nb is increased, 
we change the spatial extension of the local orbitals (Eq. 
[T5)l spanning the correlated subspace, so that the DMFT 
treatment does not apply to the same objects. Conver- 
gence studies with fixed Xm's could also be performed, 
but this is not the main scope of this article. 

Instead, we would like to emphasize some of the physi- 
cal issues when applying DMFT to different local orbitals 
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FIG. 9: (Color online) Spectral function of t^g Bloch states for 
SrVOs within the PLO scheme of LDA-hDMFT. For Nb = 12 
and Nb = 21: U=6eV. For Nb = 3: U=4eV. 



with different degrees of spatial localization. As we have 
seen, for more localized orbitals (larger TV^), we have to 
increase the value of the on-site U on the vanadium site 
to get consistent results. This is of course expected phys- 
ically. We see however that the main features of the spec- 
tral function (see Fig. [5]) are almost identical for iVfc=12 
and iVf,=21 with similar values of U. The lower Hubbard 
band (between -1.7 eV and -2.0 eV), is not far from the 
position found experimental y^^'^^i^^i^^i^^ ( between -1.5 
eV and -2.0 eV ). The Hubbard band (between 2.5 and 
3eV) is also in the range of experiments (2.5 eV). 

In the results that we have obtained, the lower Hub- 
bard band is rather systematically at an energy in which 
oxygen states already give a sizeable contribution in the 
total density of states. Experiments, however, seem to 
suggest a somewhat larger separation. We believe that 
this is due to the fact that the relative location of oxy- 
gen and y-t2g state is not accurately obtained at the 
LDA level, and that a better starting point (such as 
GW) is required to handle this problem with better ac- 
curacy. Indeed, we have verified that the intensity of the 
lower Hubbard band, and especially, of the upper Hub- 
bard band (as revealed in the Bloch-resolved t2g spec- 
tral function) is very sensitive to the precise value of the 
double-counting correction which is used. This is shown 
in Fig. llOi in which we have slightly shifted downwards 
the O— 2p states with respect to the V-i2g ones by choos- 
ing purposedly a smaller value of the double counting 
correction (3.6 eV instead of 6.6 eV). The more isolated 
the t2g states are, and the more prominent are correlation 
effects within that band (for a given U). 

In future work, calculations with extended basis sets 
should naturally face the issue of calculating the on-site 
Coulomb interaction from first-principles, but also of tak- 
ing into account Coulomb repulsion terms Upp on the 
oxygen sites, as well as inter-site repulsions Upd between 
vanadium and oxygen. We note that treating the lat- 
ter in the Hartree approximation precisely brings in a 
correction to the relative position of oxygen and vana- 
dium states, of the same nature than the double-counting 
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FIG. 10: (Color online) Spectral functions of Oxygen-p, and 
Vanadium-t29 states for U=6 eV within LDA+DMFT using 
the PLO scheme (A*';, — 12). Two values of the double count- 
ing shift are used: 6.6 eV (AMF) and 3.6 eV. 



FIG. 11: (Color online) Projected /3-NiS structure with a local 
Ni(eg) orbital obtained from a MLWF construction using the 
block of 16 Ni(3ci)/S(3p) bands. The c axis is identical to the 
vertical axis. 



terms. 

The general conclusion of this study of SrVOa at this 
stage is that our formalism is able to describe the main 
features of the experimental function within a general 
formalism which can take into account all states in the 
basis. 



IV. APPLICATION: /3-NiS 

The hexagonal form of nickel sulfide (/3-NiS) has at- 
tracted a lot of interest over the years since it exhibits 
a first-order electronic phase transition^ at about 260 
K. Whereas the high-temperature phase may be classi- 
fied as a paramagnetic metal, below the transition /3-NiS 
shows antiferromagnetic order and the resistivity behav- 
ior corresponds to characteristics of a semi metal2^ or 
a degenerate semiconductor—. Hence the term "anti- 
ferromagnetic nonmetal" is commonly used for the low- 
temperature phase. Note that hexagonal NiS is only 
metastable at room temperature, the true stable phase 
is given by the millerite structure-^. The crystal struc- 
ture (Fig. [11]) of /3-NiS is of the NiAs-type (space group 
Pd^/mmc) with two unit cells in the primitive cell^. 
In this rather simple structure the NiSg octahedra share 
edges within the ab plane and share faces along the c 
axis. There is a slight decrease of the cell parameters 
((5a/a~0.3%,(5c/c~l%) below the transition, giving rise 
to a volume collapse of the order of 2 percent. Our aim 
in this work is not to perform a detailed study of the 
metal-to-nonmetal transition of /3-NiS, since this would 
involve deeper considerations concerning the role of mag- 
netism above and below the transition. For an overview 
of this, still controversially discussed, topic see e.g. the 
review by Imada et ali^. 

Here we mainly want to use the high-temperature 
phase of /3-NiS as an example for a correlated metal with 
strongly hybridized Ni(3d) and S(3p) states, formally in 
the charge-transfer regime of the Zaanen-Sawatzky- Allen 



classification scheme^. In this regard it will become clear 
that this compound may not easily be treated in the "tra- 
ditional" LDA+DMFT scheme via projecting solely onto 
low-energy states close to the Fermi level. Thus the goal 
is to use the here outlined projection technique of in- 
terfacing LDA with DMFT in order to explore the im- 
portance of electronic correlations for the local spectral 
function. 



A. LDA investigation 

Electronic structure calculations for /3-NiS date back 
to the original work of Mattheis^. Here the DFT-LDA 
calculations are performed with a mixed-basis pseudopo- 
tential (MBPP) code2i described in appendix lA 21 For 
the lattice parameters we use a= 3.440 A and c=5.351 
A from Ref. IM Fi gure [T2] shows the resulting band 
structure and density of states (DOS) from this compu- 
tation. It is seen that the common block of hybridized 
Ni(3(i) and S(3p) bands is isolated from the remaining 
bands. The dominantly S(3s)-like bands are ~5.5 eV 
below in energy^^-, while the bands above starting with 
mainly Ni(4s) character are separated by ~1.2 eV from 
that block. The latter has a width of ^--^9 eV and consists 
of three prominent peaks below the Fermi level. The 
dominant d^L peak^ at -1.5 eV (Pi), a second at -4.4 
eV (P2) and the last at -6.5 eV (P3). Whereas P2 has 
mixed Ni(3d)/S(3p) character, P3 is a nearly pure S(3p) 
peak and does not stem from bonding between Ni and S. 

In order to obtain more detailed information about the 
involvement of the different orbital sectors, the orbital- 
resolved DOS is additionally incorporated in Fig. [T^ . 
Furthermore, the fatband representation of the decompo- 
sition of the Bloch bands is presented in Fig.[T31 Because 
of the hexagonal symmetry the Ni(3d) multiplet splits 
into two degenerate Cg levels, two degenerate Eg levels, 
and one Aig level per atom, respectively. As usual, the 
Cg states hybridize more strongly with the S(3p) states 
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than the remaining {Eg,Aig) states. Due to the increas- 
ingly filled d shell of Ni, the {Eg,Aig) levels are more- 
over also nearly completely occupied in /3-NiS. In num- 
bers, the LDA filling, including degeneracy, amounts to 
(2.9;3.9,1.9) for {eg;Eg,Aig). The fatband representation 
clearly shows that the hybridization between the differ- 
ent orbitals are strong and its not quite obvious to single 
out a low-energy regime for this compound within the 
central block of bands. Although the bands at the Fermi 
level appear to be dominantly of Cg character, the corre- 
spondig orbitals have strong weight at lower energy, too. 
Despite the strong filling of the {Eg,Aig) states, abandon- 
ing those orbitals in a minimal hamiltonian description of 
/3-NiS is likely to fail. Fluctuations originating from hole- 
like states may be important in the end to reach a bet- 
ter understanding of the complex magnetic behavior of 
this material. Note also in this respect that the effective 
bandwidth of the energy range where the eg character is 
manifest is roughly twice as large as the corresponding 
energy range for {Eg,Aig). Hence this reduced relative 
bandwidth for {Eg,Aig) within the Ni(3d) multiplet will 
also have influence on the degree of the orbital-resolved 
correlation effects. 
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FIG. 12: (Color online) LDA data for /3-NiS. (a) Band struc- 
ture, (b) DOS. For the local Ni(3d)/S(3p)-DOS the cutoff 
radius was half the nearest-neighbor distance, respectively. 
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FIG. 13: (Color online) LDA fatband decomposition for /3- 
NiS. (red/gray) m{Eg) and (blue/dark) Ni(^i3). The local 
Ni(3d)/S(3p) projection cutoff radius was half the nearest- 
neighbor distance, respectively. 



B. LDA+DMFT investigation 

The influence of correlation effects in /3-NiS has been 
investigated, experimentally and theoretically, by several 
author o^°i'^-^''*^i'^^i^'^i'^^ . It is generally believed that this 
compound is moderately correlated in view of the strong 
hybridization of the ~3/4 filled eg states with sulfur. A 
value for the Hubbard U of the order of 4—5 eV was es- 
timated from modelling the measured spectral function 
derived from photoemission experiments^i^i^. These 
experiments in the metallic phase are in rough agreement 
with the shown LDA DOS below the Fermi level in so far 
as they also reveal three peaks with comparable relative 
intensity as the theoretical set (Pi, P2, P3). However, an 
effective single-particle picture appears to be insufficient 
to understand those peaks, especially when going to the 
nonmetallic phased. It is theoretically expected that a 
lower Hubbard band, i.e., satellite, originating from the 
Ni(3(i) states should be located within the Bloch states 
deep in energy (starting around 6 eV) and dominantly 
characterized by S(3p). This idea relies on the fact that 
/3-NiS may be viewed as belonging to the charge-transfer 
category of transition-metal chalcogenides, yet not be- 
ing that strongly correlated for the lower Hubbard band 
appearing completely below the S(3p) states. 

The strong hybridization between Ni(3d) and S{3p) 
renders the PLO version of the LDA+DMFT method 
most suitable for this compound. Concerning the value 
of U we chose a pragmatic approach and performed the 
calculations for two possibly reasonable choices, i.e., ?7=4 
eV and U =5 eV. The value of J is certainly less materials 
dependent and was fixed to J=0.7 eV. In order to take 
care of the double counting, we used the formalism of 
fixing the local total charge (see appendix [B]) . Note that 
for the present crystal structure there are two symmetry- 
equivalent Ni atoms in the primitive unit cell, i.e., two 
correlated sites R. The local Green's function and self- 
energy were thus computed by symmetrizing the site- 
resolved quantities. For the projection onto local orbitals 
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we limited the number of bands to iVb=16, i.e., all bands 
of the central block around the Fermi level are used. 
This renders the corresponding correlated subspace al- 
ready rather localized, e.g., the effective Ni(3(i)-like WFs 
from this set of bands are not expected to leak much to 
the sulfur sites (see also Fig. [TT|). 

By explicitly including the correlation effects, the 
orbital- resolved fillings in the 3d shell of Ni do not change 
relevantly, thus effects due to changes in orbital popula- 
tions induced by correlations are not expected to play an 
important role for this compound. Figure [TH exhibits the 
resulting local spectral function for the different values of 
U at inverse temperature (3=10 eV~^. It is seen that the 
influence of correlation effects on the metallic spectral 
function are indeed rather subtle for this compound. A 
lower Hubbard band appears to show up for C/>4 eV 
within the dominant S(3p) energy regime. It is however 



hard to extract this atomic-like excitation from the to- 
tal spectral function (as in experiment'"'"'*^). The plot 
of the impurity spectral function in Fig. [15] reveals that 
the {Eg,Aig) orbitals yield indeed effective bands with 
smaller bandwidth than the eg orbitals. It seems that 
the contribution to the lower Hubbard band then also 
stems more significantly from the former set of orbitals. 

In the end, the total spectral function computed for 
C/=5 eV shows close resemblance to recent experimen- 
tal curves obtained from photoemission^. Note however 
that it appears tricky to disentangle the lower Hubbard 
band from the P3 peak of in principle pure S(3p) content. 
Our calculation suggests that the correlation effects have 
also impact on the lowest S(3p)-like states. Further stud- 
ies are necessary to clarify this issue. 



V. CONCLUSIONS 
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FIG. 14: (Color online) The LDA-I-DMFT local spectral func- 
tion for /3-NiS derived from the local Green's function in Bloch 
basis at (3=10 eV~^. The contribution from the upper 10 
bands of the Ni(3d)/S(3p) block was encoded red, the one 
from the lower 6 bands was encoded blue. This guidance to 
the eye should roughly seperate dominant Ni(3d) from domi- 
nant S(3p) character. 
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FIG. 15: (Color online) The LDA-f-DMFT local impurity 
spectral function for /3-NiS at /3=15 eV^^. 



We have implemented an effective and flexible 
LDA-I-DMFT scheme into two electronic structure for- 
malisms based on a plane-wave description, i.e., pro- 
jector augmented- wave and mixed-basis pseudopotential. 
The orbitals defining the correlated subspace C, in which 
many-body effects are included, are constructed by pro- 
jecting local atomic-like orbitals onto a restricted set W 
of Kohn-Sham Bloch states, similar to the procedure 
adopted in Ref. h Orthonormalization of the projected 
orbitals yields effective Wannier functions spanning C. 
These WFs are not unique: they depend on the energy 
window covered by the Bloch functions in W (the larger 
the window, the more localized the WFs are). Although 
more sophisticated constructions^iii of explicit WFs will 
surely remain an important tool in the LDA-I-DMFT con- 
text, we feel that the more straightforward projection 
technique will render future developments in this area 
easier and more efficient. 

Using this method, we have investigated SrVOa and 
compared different implementations involving different 
choices for the set of Bloch states W and for the cor- 
responding local orbitals spanning C. We have shown 
that the basic physical findings for this compound are 
consistent with previous LDA-I-DMFT treatments, inde- 
pendently of the chosen implementation. However, the 
present study allows for an explicit treatment of ligand 
states, which raises several issues which should be the 
subject of further studies. One of these issues is a first- 
principle determination of the on-site Coulomb matrix el- 
ements in the different choices of local orbitals and basis- 
set. Our study supports the expected fact that a larger 
value of Udd has to be taken when a larger energy range 
(and more localized orbitals) is considered. Furthermore, 
a proper treatment of the on-site repulsion on oxygen 
sites Upp, as well as of oxygen-transition metal inter-site 
repulsion Upd should certainly be considered. This is- 
sue is tightly connected to the choice of double-counting 
correction, which, as we have shown, plays a significant 
role. 
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We have also presented a first LDA+DMFT study of 
metallic /3— NiS, a charge-transfer compound for which 
the inclusion of ligand states is crucial. More detailed 
studies, extending also into the nonmetallic regime, are 
needed in order to gain more insight into this "traditional 
companion" of the more famous NiO compound where 
such investigations recently took place^i^. 

In general, the explicit inclusion of the ligand states in 
the LDA+DMFT description, based on generic highly- 
accurate electronic structure codes, will open the door 
to new possibilities for the investigation of strong- 
correlation effects in real materials. 



APPENDIX A: CALCULATION OF THE 
PROJECTION P^^ 

In this appendix, we describe the details of the calcula- 
tion of the projection P^^^ from KS Bloch states onto lo- 
cal orbitals in both the PAWi^ and the MBPPi^ method. 

1. Projected Augmented- Wave (PAW) method 

In the projected augmcnted-wave formalism, two kinds 
of atomic functions are used: the pseudowavefunction ip 
and the true wavefunction ip. They are used to recover 
the correct nodal structure of the wavefunction near the 
nucleus (see Eq. (9) of Blochl et a^). A direct cal- 
culation of the projection P;^j^(k) would thus require to 
compute the three terms resulting from this equation, 
namely (pi is the projector for angular momenta h, 
and its projection rrii): 

^'m.(k) = (X^l^k.) +^(K|$k.)((Xmk.> - 

(Al) 

However, atomic d or f wavefunctions arc mainly local- 
ized inside spheres. As a result, we can compute the pro- 
jection only inside sphere, in the spirit of current LDA+U 
implementation^!^ in PAW. In this particular case, the 
first and third terms of Equation lAU cancel with each 
others. This cancellation is only exact for a complete set 
of projectors in the energy range of the calculation. This 
completeness can be easily tested during the construc- 
tion of projectors and partial waves. The projection thus 
writes as a sum over projectors p„^: 

P;^,(k)=^(p„J$k.)(x^|^n,) (A2) 

We note that integrals in this equation are computed 
only inside sphere. It implies that the projection is done 
on an unnormalized wavefunction %. Nevertheless, the 
later normalization of the projection P,^j,(k) will make a 
normalization of x redundant. This implementation has 
been made with the code ABINIT^i^^. 

Additionally, we implemented the proposed 
LDA+DMFT scheme in a somehow simplified interfac- 
ing with the PAW-based Vienna Ab Initio Simulation 



Package^^ (VASP) by using only the PAW projectors 
from the standard output for P;^^(k) . Although 
approximate, this latter identification yielded (after, of 
course, proper normalization) very similar results in 
comparison with the more well-defined explicit coding 
described above. 



2. Mixed-Basis pseudopotential (MBPP) method 

The combination of normconserving pseudopotentials 
with a mixed basis of plane waves and localized orbitals 
in order to represent the pseudo crystal wave function is 
the main ingredient^ of our employed mixed-basis pseu- 
dopotential code^. More concretely, the MBPP band 
structure code uses the following representation for the 
KS pseudo wave functions: 

l*k.)=E'/^kJk + G)+E/?kLlO . (A3) 
with 

(r|k + G) ^ ^e*(k+G)T (A4) 

T 

In these eqns., 7 denotes the atom in the unit cell, flc the 
volume of the unit cell, v the band and l,m are angular 
momentum and azimutal quantum number. The plane 
waves in this basis extend up to a chosen energy cutoff 
i?P™. The analytical form of the local orbitals (pjji reads 

^7/,{v')=i'f.,iKi„,{r') , r' = r-T-R^ , (A6) 

whereby AT;™ descibes a cubic harmonic and the radial 
function f^i is an atomic pseudo wave function modified 
with a proper decay and cutoff function in order for f^i to 

vanish beyond some chose radial cutoff ri''''-' . Note that 
the local orbitals are orthonormal and by definition do 
not overlap between neighboring sites, i.e., 

/ dr\f^i{r)\' = 1 (A7) 
Jo 

{4>t!:\(t>t) = S-^'^SviS^.^ . (A8) 

By identifying the local orbitals (/)^ as the objects to 
project onto (i.e. the Xm in the outlined LDA+DMFT 
formalism), the projection P^j^(k) may be written with 
the help of ((X3|) and ((X8|l as 

(A9) 
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APPENDIX B: DOUBLE COUNTING 

The double-counting shift (used in eq. [5]) is neces- 
sary in the formahsm in order to correct the fact that 
electronic correlations are already treated in an aver- 
age way within LDA. Since there is no unique way to 
extract this double counting from LDA, two alternative 
formulations were used in this work. First the so-called 
Around Mean Field^i^ (AMP) method, which is usu- 
ally adequate for metals. Second, a different formalism 
in which the double counting is such that the electronic 
charge computed from the local noninteracting Green's 
function and the one computed from the interacting im- 
purity Green's function is constraint to be identical: 

TrGi„,p(zc^„) = TrG|°p(*^n) ■ (Bl) 

Here G^^^^ is generally given by 

k cta' 

X |[*c^„ + ^-HKs(k)]-'| . (B2) 

Alternatively one may want to use the Weiss field Q{iu!n) 
instead of the local noninteracting Green's function in 
eq. (|B1[) . Note that in any case this formalism asks for 
an additional convergence parameter Sfidc to fulfill (|Bip . 
Hence this parameter has to be found iteratively together 
with the total chemical potential /x in the DMFT cycle. 

For SrV03 (with Nt = 12), the respective double- 
counting shifts are 6.6 eV (AMF) and 5.5 eV (fixing local 
charge) . However these two calculations give very similar 
results. 



APPENDIX C: MAXIMALLY LOCALIZED 
WANNIER FUNCTIONS IN PAW 

To define maximally-localizedii WFs, we need the fol- 
lowing quantity: 

Mj^:::^ = (^„,kie-'^-nvi/„,k+b) (ci) 

In the PAW framework, this quantity can be expressed 
as a function of the pscudo wave function the projec- 
tors p, the atomic wave function cj) and the pseudo atomic 



wave function (f>. We use the expression for an operator 
A wfthin the PAW formahsm (Eq. (11) of Blochl et a!^^). 
We thus obtain: 

+ ^{^Km\Pi){P3\^k+b,m) (C2) 

X ((0,|e-*-"-|</',)-(0.|e-''-'-|0,)) 

A similar expression has been used for Ultra-Soft Pscu- 
dopotentials by Ferretti et a^. The first term is com- 
puted in the Fourier basis, whereas the second is com- 
puted in the radial grid. We use the expansion of e"'*' '" 
in spherical harmonics and bessel functions in order to 
compute its expectation value over atomic wavefunctions. 

The overlap matrix Mmju"^ is used to minimize the 
spread and localize the WFs. This is done with the pub- 
licly available Wannier90 code^. 

The code is tested on t2g MLWF for SrV03 using only 
the three KS t2g bands. The spread of one of the WFs 
for a 8 X 8 X 8 k-point grid centered on the P point with 
our implementation is 6.96 a.u^. In excellent agreement 
with the FLAPW result (6.96 a.u^ in Ref.H). 
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Note added: As the writing of this work was being 
completed, we became aware of the related work of Ko- 
rotin et al. (larXiv:0 801.3500) also reporting on an imple- 
mentation of LDA-I-DMFT in a pseudopotcntial plane- 
wave framework. 
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